Condensation in Temporally Correlated Zero-Range Dynamics 
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Condensation phenomena in non-equilibrium systems have been modeled by the zero-range pro- 
cess, which is a model of particles hopping between boxes with Markovian dynamics. In many 
cases, memory effects in the dynamics cannot be neglected. In an attempt to understand the pos- 
sible impact of temporal correlations on the condensate, we introduce and study a process with 
non-Markovian zero-range dynamics. We find that memory effects have significant impact on the 
condensation scenario. Specifically, two main results are found: (1) In mean-field dynamics, the 
steady state corresponds to that of a Markovian ZRP, but with modified hopping rates which can 
affect condensation, and (2) for nearest-neighbor hopping in one dimension, the condensate occupies 
two adjacent sites on the lattice and drifts with a finite velocity. The validity of these results in a 
more general context is discussed. 
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The Zero-Range Process (ZRP) is a paradigmatic 
model for mass transfer in non-equilibrium systems [l], Q . 
In this model, particles hop between boxes with hopping 
rates u(n) which depend only on the occupation number 
n of the departing box. The steady state distribution 
function of the box occupation numbers factorizes into 
a product of single-box terms, which makes this model 
amenable to theoretical studies. For rates which for large 
n are of the form 



u(n) = 7(1 + b/n), 



(1) 



it has been shown that the model exhibits a condensa- 
tion transition at high densities as long as b > 2. At 
densities above the condensation transition, one of the 
boxes is occupied by a macroscopic number of particles, 
constituting a real-space condensate. The parameter 7 
sets the time scale of the process. 

The ZRP has been used to study condensation in 
shaken granular gases where the shaken particles can 
move between compartments. Under certain conditions 
condensation is observed, whereby most particles accu- 
mulate in one compartment [H, Q] • Similar condensation 
phenomena may occur in network rewiring, resulting in 
a hub with a macroscopic linking number [B|, Q. The 
ZRP has also been used to model phase separation in 
driven diffusive systems and traffic flow models, where, 
in the steady state, high and low density regions coexist 
0, H H E3] • In modelling one dimensional driven diffu- 
sive systems, the length of, say, a high density region in 
the driven system corresponds to the occupation number 
of a box in the ZRP. The flow of particles between high 
density domains corresponds to particle hopping in the 
ZRP. 

In reality, mass transfer processes usually constitute a 
coarse-grained description of more complex microscopic 
dynamics. This frequently results in temporally corre- 
lated dynamical processes, which are not captured in the 
Markovian dynamics of the ZRP ll|. It would thus be 



on the collective behavior of the ZRP, particularly on the 
occurrence of condensation. 

In this Letter, we consider a non-Markovian ZRP, 
where the hopping processes depend not only on the 
present occupation of a site but also on jump events that 
have occurred in the past. We find that this has sig- 
nificant impact on condensation. Specifically, two main 
results are derived. The first is that memory effects on 
steady state condensation can be captured by a Marko- 
vian ZRP with an effective hopping parameter 6 c ff ^ b 
which controls the condensation transition. The second 
effect is that in the case of asymmetric nearest-neighbor 
hopping on a ring of L sites, the condensate drifts with 
a finite velocity which scales as 1/L. 

The non-Markovian ZRP introduced in this work con- 
sists of L boxes i, each characterized by the particle num- 
ber rii and an internal "clock" Tj. The total number of 
particles is N . The clock proceeds irregularly in integer 
steps and is reset to zero each time a particle jumps onto 
site i, thus keeping a memory of the history of the pro- 
cess. A configuration of the system is then given by the 
set (n, t) = {(ni,Ti)}f =1 . The hopping rates out of site 
i are taken to depend both on the occupation of the site 
rii and on the state of the clocks Ti, but not on the occu- 
pation of any other site, in accordance with the general 
approach of the ZRP. 

Here we consider the following dynamics: particles hop 
between sites with rate u(n, t). The internal clock on the 
target site is reset to zero together with a jump. Inde- 
pendently of the jump processes, at each site the internal 
clock is incremented by one unit with a constant rate 
c. For a jump from i to j, these two processes can be 
schematically summarized by 



{ni,Ti), (nj,Tj) > (ni-l,Ti), {nj + l,Tj = 0) 



(nun 



-> (m,Ti+i), 



(2) 



of interest to explore the effect of temporal correlations 



The particle jump process, when taken by itself, is non- 
Markovian, since the rate of a jump depends on how 



2 



much time has past since a particle last hopped into 
the jump site. Choosing the target site j ^ i uniformly 
(i.e., considering a fully connected graph) corresponds 
to mean- field (MF) dynamics, while restricting the tar- 
get site to j = i + 1 corresponds to totally-asymmetric 
nearest-neighbor hopping on a ring. 

Notice that the full process of particle jumps and 
clock increments together defines a Markovian process. 
Therefore, the model may be implemented by a discrete- 
time Monte-Carlo version of these dynamics with ran- 
dom sequential update which is defined as follows: let 
Pmax = max niT [u(n, r) + c\. For the Monte-Carlo up- 
date pick a pair of sites uniformly and attempt to 
make one of the following changes: (i) move a particle to 
the target site j with probability u(rii, Ti)/p max and reset 
the clock on the target site j to zero, or, (ii) increment 
the internal clock on the starting site i with probability 
c/pmax- A total of L consecutive updates constitute one 
Monte-Carlo sweep. 

Let us first consider the MF dynamics, whereby the 
target site j is chosen uniformly. With this dynamics 
the steady state distribution factorizes into a product of 
single site terms in the thermodynamic limit. The single- 
site occupation and clock probability in the steady state, 
P(n, t), can then be found by examining a single site with 
a "mean-field" incoming current J generated by all other 
sites. This is equivalent to a solution of a single site in 
the grand-canonical ensemble. The master equation for 
the single site probability is 



dP{n,r) 
dt 



-P(n,r)[J + c + u(n,T)] + 



JP(n-l)S Ti0 + cP(n,r-l) + u(n+l, r)P(n+l, r), (3) 
where the single site marginal distribution is defined by 

P(n) EE J2 T P ( n i T )> and J = J2n.T U ( n ' T ) P ( n ' T ) is tlic 

incoming current. These equations are valid also for n — 
and r = by defining P(-l,r) = P(n, -1) = and 
w(0,t) = 0. 

In the steady state, where dP(n 7 r)/dt = 0, one obtains 
by summing Eqs. Q over r 

u{n)P(n) = JP(n-l). (4) 

Here u(n) is the mean hopping rate out of a site with n 
particles 

_ Er p (n,T)u(n,T) 



u(n) 



(5) 



Eq. f4| expresses the balance between the probabilities 
to hop into and out of a site with n particles. The 
marginal distribution P(n) is therefore the same as that 
of a Markovian ZRP [12j |. but with an effective hopping 
rate u(n) . This gives the steady-state occupation proba- 
bility 



P(n) = P(0)J n f(n) with f(n) 



i=l 



(6) 



Here P(0) _1 = 1 + Y^Li J n f( n ) ensures proper nor- 
malization of P(n). The occurrence of condensation is 
determined by the asymptotic behavior of P(n). For 
P(n) ~ n~ a condensation takes place for a > 2 0]. In 
the Markovian case, ones has a — b for hopping rates of 
the form |T]). 

To analyze condensation in the non-Markovian case 
we proceed by considering, for simplicity, hopping rates 
of the form 



u(n, t) 



(7) 



r = ("off" state) 
u(n) t > 1 ("on" state). 

In this case, whenever a particle hops into a site the site 
is switched to an "off" state in which no particles can hop 
out. Only after its clock reaches r = 1 the site is turned 
"on" again, in which case particles hop out with a rate 
u(n). This special case will be called the on- off model. 
In this case, the clock has in effect only two states: t = 
corresponding to "off" , and r > 1 corresponding to "on" . 
Correspondingly, the state of the site can be character- 
ized by P s{n) = P(n, 0) and P on (n) = J2 T >i p ( n > T )- 

The master equation for the stationary probability dis- 
tribution ([3]) is then 

P oS (n)[J + c] = JP(n-l) (8a) 
Pon(n)[J + u(n)] = cP oS {n) + P on {n + l)u(n+l). (8b) 

First we study the probability P g to find a site in the 
"off" state. By summing over n, Eq. ([8a| yields 

J 



Pn 



C+J' 



Further, from Eqs. ©, (5aJ) and © we find 

PoffH 



u{n)P(n) P off _ 

— T- = —;-u{n)P(n). 

J + c J 



(9) 



(10) 



In addition, for u(n, t) of the form 0, Eq. gives 

u(n) 



PonH 



u(n) 



P(n). 



(11) 



As discussed above, the occurrence of condensation de- 
pends only on the asymptotic behavior of the effective 
hopping rate, u(n). The effective hopping rate follows 
from ([TO]) and (TTTj) using P e(n) + P on {n) = P(n), and is 
given by 

u[n) J uyn) 

This states that the mean time between hops from a site 
with n particles is equal to the mean time this site is in 
an "off" state plus the time it takes a particle to hop out 
once the system is already "on" . 

In the case of u(n) of the form ([1]), it can be seen from 
(fT2]) that to leading order in 1/n the effective hopping 
rates are are again of the form (fTJ) , given by 



u(n) 



J 



IV n 



(13) 
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FIG. 1: Snapshots of the on-off model with totally- 
asymmetric nearest-neighbor hopping on a ring, showing the 
occupation numbers m at and in the vicinity of the conden- 
sate at four points in time. Here, L = 1000 sites, p = 10 
b — 5.5 and c = 1. The condensate occupies two sites and 
drifts with a constant mean velocity. 
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FIG. 2: The position of the most occupied site i max and its 
occupation number n max as a function of time, on a ring of size 
L = 1000 with totally-asymmetric nearest-neighbor hopping. 
Simulation parameters: p — 10 b = 5.5 and c = 1. 



with 



7 c + J L L 
o c ff = ; r < 6. 



c+ J+ 1 
From ©, we then have 

"J(J + c+ 1) 



P(n) 



J + c 



"beff 



(14) 



(15) 



To analyze the condensation transition we note that 
for J(J + c + 1)/(J + c) < 1, the distribution decays 
exponentially with n and the system is in a subcriti- 
cal homogeneous phase. For a critical current given by 
J c = (c + J c )/(c + J c + 1), the decay is algebraic, and 
one has condensation for b c ff > 2. For the critical cur- 
rent we find 



2 



1 + --1 

c 



which allows us to write 



b c g = J c -b. 



(16) 



(17) 



Therefore, condensation takes place for hopping param- 
eter b > | + 4/c — l^j which is larger than 2, in 
contrast with the Markovian case for which the critical 
value for condensation is b = 2. This means that for 
densities p = N/L above a critical density p c , the av- 
erage number of particles in the condensate is Ar con d = 
L(p — p c ). All other sites have an average "background" 
density p c . 

We now turn to a ring geometry with nearest neighbor 
hopping and fully asymmetric dynamics. In this case, 
the stationary distribution does not factorize. In order 
to get insight into the condensation mechanism in this 
case we carried out numerical simulations of the on-off 
model with u(n) given by Eq. (JTJ) . 



Simulations in the condensed phase indicate that the 
condensate drifts with a finite velocity. In addition, we 
find that in contrast to previously known condensation 
phenomena, the condensate typically occupies two adja- 
cent sites i and i + 1. These observations can be ex- 
plained by a more detailed microscopic investigation of 
the dynamics. It turns out that the drift takes place via 
a "slinky" motion where particles hop from site i to i + 1 
at a rate u(rii) which is approximately constant at large 
rij, leaving site i + 1 in a predominantly off-state. Thus 
particles accumulate at site i + 1 until site i is emptied, 
giving the clock at i + 1 the chance to reach the on-state 
for durations of time sufficiently long to allow particles 
to escape. Then particles start to hop from site i + 1 to 
site i + 2 in the same fashion. 

This mechanism suggests that the drift velocity w CO nd 
is inversely proportional to the number of particles in the 
condensate N CO nd, i- e -i 



N-N c = L(p-p c ). 



(18) 



In the thermodynamic limit, the velocity of the conden- 
sate vanishes. Superimposed on this motion, the con- 
densate can melt and reappear at some other site of the 
lattice, similar to what happens in the Markovian case. 
This happens on a characteristic time which scales with 
the system size to a power larger than 2 13|, [l4|, [15 1 . 

In Fig. [T] we present snapshots at different times of the 
on-off model with totally asymmetric nearest neighbor 
hopping for L = 1000. One clearly sees that the con- 
densate occupies two adjacent sites with varying relative 
occupation of the two sites, consistent with the slinky 
motion described above. The drift of the condensate is 
evident in the figure. In order to demonstrate the slinky 
motion in more detail, we present in Fig.[2]a plot showing 
the position of the most occupied site i max and its occu- 
pation number, n max , as a function of time. The occupa- 
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FIG. 3: The occupation probability P(n) of a single-site in 
a ring with totally-asymmetric nearest-neighbor hopping and 
different densities, as obtained from Monte-Carlo simulations. 
Here L = 1000, and curves for subcritical density (p = 3), 
supercritical density (p = 10) and at the critical region (p — 
4.1) are presented. The horizontal line indicates P p i a toau of 
(|19|) . where p c was obtained from the simulation. The inset 
shows a similar plot of P(n) for a Markovian ZRP of L = 1000 
sites with b — 5, in the subcritical (p = 0.5) critical (p — 1) 
and supercritical (p — 4) phases. 



tion number n max oscillates in time with approximately 
constant frequency. Typically it decreases linearly until 
it reaches its minimal value, when i max increases by 1 
and n max starts increasing. 

In Fig. [3] we present the single-site occupation prob- 
ability distribution P(n) for various densities. At high 
densities the distribution exhibits a plateau which reflects 
the particle distribution among the two sites which con- 
stitute the condensate. This is in contrast with a Marko- 
vian ZRP where the condensate is supported by a single 
site, which results in a sharp peak in P(n). The value 
of P(n) at the plateau in the non-Markovian case may 
be estimated for p above the critical density and large L 
using the slinky motion of the condensate. The probabil- 
ity that a given site carries the condensate is 2/L, and in 
such a site there is an approximately uniform probability 
to find any occupation < n < iV con d = L(p — p c ). Thus, 



P plateau 



2 
L 



1 



Hp- Pc)' 



(19) 



This estimate agrees well with the plateau value in Fig. [3] 
For small densities, P(n) decays exponentially, indicating 
the absence of a condensate. For the systems size studied 
in this figure, the distribution at small values of n does 
not allow to extract a power law decay as expected for 
the condensation transition. At density p = 4.1 there 
is a range of n for which P(n) seems to follow a power 
law with b e g ~ 4. This value differs significantly from 
b = 5.5, which is the expected power for the Markovian 
ZRP. 

So far, we discussed in detail a simple on-off model 



which has been used to demonstrate the effect of non- 
Markovian dynamics on the steady state distribution 
and on condensation. For MF dynamics, the asymp- 
totic behavior of the effective rates can be obtained for 
the more general class of models with rates of the form 
u(n,T) = u(ji)v(t). Generally, 6 c ff may be larger or 
smaller than the "bare" value of b. Interestingly, one can 
introduce an on-off model with nearest-neighbor hopping 
for which the stationary distribution factorizes strictly 
even on finite lattices and where the effective hopping 
rates can be computed exactly. In this model particles 
jump with rates (JT]) but the advancement of the clock 
depends on the clock states of neighboring sites [l6| ■ 

The findings of this work suggest that the tempo- 
ral correlations may significantly alter the condensation 
transition and the nature of the condensate in general 
driven systems. It would be of interest to study broader 
classes of temporal correlations in order to explore other 
scenarios of non-equilibrium condensation. 
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